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POLARIZED LINE PROFILES AS DIAGNOSTICS OF ClRCUMSTELLAR 
GEOMETRY IN TYPE UN SUPERNOVAE 
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RESUMEN 

Favor de proporcionar un resumen en espanol. If you cannot provide a Spanish abstract, the 
editors will do this. Supernovae of type Iln possess spectral signatures that indicate an intense interaction 
. between the supernova ejecta and surrounding dense circumstellar material cast off by the star in pre-explosion 
mass-loss episodes. Studying this interaction can yield clues to the nature of Type Iln progenitors and their 
mass loss history. In particular, polarization spectra of Type Iln's show complex line polarization and position 
angle features that arise from a combination of geometrical and optical effects. I have constructed a Monte 
O , Carlo code that simulates the transfer of the Ha line through circumstellar shells with various geometrical 

Q configurations and optical characteristics. The superposition of broad and narrow line components produced 
in different regions of the circumstellar environment and modified by electron and line scattering, hydrogen 
absorption, thermal emission, and geometrical and viewing angle effects gives rise to a variety of polarized 
line shapes in the model spectra. Comparison of these results with recent high-quality spectropolarimetric 
observations of Type Iln supernovae suggests that a model "shock" region between the supernova photosphere 
and the circumstellar shell is necessary to produce the narrow polarized emission features at the rest wavelength 
of Ha seen in some Iln's. Further model results point toward other features in the polarized line profile that 
■<sj- , can be used to constrain the characteristics of the circumstellar material in these intriguing objects. The code's 
usefulness will be extended by the treatment of Dopplcr effects due to expansion of the circumstellar scattering 
region, such as those that characterize the polarized Ha profiles of the Type Iln SN 1997eg. 
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ABSTRACT 

Supernovae of type Iln possess spectral signatures that indicate an intense interaction between the supernova 
ejecta and surrounding dense circumstellar material cast off by the star in pre-explosion mass-loss episodes. 
Studying this interaction can yield clues to the nature of Type Iln progenitors and their mass loss history. 
In particular, polarization spectra of Type Iln's show complex line polarization and position angle features 
that arise from a combination of geometrical and optical effects. I have constructed a Monte Carlo code that 
simulates the transfer of the Ha line through circumstellar shells with various geometrical configurations and 
optical characteristics. The superposition of broad and narrow line components produced in different regions 
of the circumstellar environment and modified by electron and line scattering, hydrogen absorption, thermal 
' emission, and geometrical and viewing angle effects gives rise to a variety of polarized line shapes in the model 
spectra. Comparison of these results with recent high-quality spectropolarimetric observations of Type Iln 
supernovae suggests that a model "shock" region between the supernova photosphere and the circumstellar 
shell is necessary to produce the narrow polarized emission features at the rest wavelength of Ha seen in 
some Iln's. Further model results point toward other features in the polarized line profile that can be used 
to constrain the characteristics of the circumstellar material in these intriguing objects. The code's usefulness 
will be extended by the treatment of Doppler effects due to expansion of the circumstellar scattering region, 
such as those that characterize the polarized Ha profiles of the Type Iln SN 1997eg. 
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DIATIVE TRANSFER — STARS: MASS LOSS — SUPERNOVAE 

1. INTRODUCTION collapse types of these stellar e xplosions are inher- 

ently aspherical phe nomena (e.g., Kozma et al.ll2005t 
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Burrows et al.ll2006l ) . Observations of net continuum 



Our understanding of supernovae (SNe) has 

broadened in recent years to include the recogni- 
, . ,i , , , T ,i i i ,i polarization m both supernova (SN) types have pro- 

tion that both the thermonuclear and the core- . , , , . , r . . . ' * ri1kT 

vided key e vidence of the intrinsic asymmetr y ol SN 
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SNe also display line polarization features in ad- 
dition to broadband continuum polarization; these 
line effects are often more complex than simple de- 
polarization by complete scattering redistribution, 
and they can provide specific clues to the nature of 
SN ejecta and their surrounding circumste llar me- 
dia (e.g., [Wang| |200i EeonaxdeEal] Ex)!). How- 
ever, line polarization can be produced by a combi- 
nation of many different optical and geometrical ef- 
fects, so its interpretation is not straight forward. De 



tailed radiative transfer modeling (as in lKasen et al 



20031 ) is often the best way to understand the ef- 
fects that give rise to polarized lines in SN spectra. 
In th e case of Type Iln ( "narrow- line" ; iFilippenko 
19971 ) supernovae, the situation is complicated by 
the presence of circumstellar material (CSM) sur- 
rounding the SN ejecta that becomes excited by the 
UV and X-ray photons from the SN explosion. Line 
polarization signatures in Type Iln SNe are super- 
positions of those arising from the ejecta and those 
arising from the CSM. Disentangling the two can 
be difficult but worthwhile, as it allows us to probe 
the nature and structure of the CSM to an extent 



not possible with sp ectroscopy alone ([Leonard et al 



not possible witn sp ' 

2000al IWandl200ll ). Since the CSM of a Type Iln 
supernova most likely represents material ejected by 
the progenitor star, studying its geometrical and op- 
tical characteristics provides a link to the mass-loss 
episodes and stellar winds of massive stars in their 
late stages of evolution. 

2. CODE 

I have developed a Monte Carlo radiative transfer 
code called SLIP (for Supernova Line Polarization) 
that simulates the ways polarized line profiles are 
created in Type Iln supernovae. SLIP uses the 
three-dime nsional spherical p o lar gr id structure de- 
scribed bv IWhitnev fc Wolfj (|2002l ); it tracks vir- 
tual photons as they arise from a model SN photo- 
sphere and scatter in a circumstellar density distri- 
bution with wavelength-dependent emission, absorp- 
tion, and scattering characteristics. Similar codes 
have enjoyed success in analyzing related scenarios 
such as h ot star envelop es with aspherical wind ge- 
ometries (IHarriesI 120001 ) and ejecta-hole c onfigura- 



tions in SNe Type la (jKasen et all 120041 ). Unlike 



most previous cod es that treat SN line polarization 
re.g.. lHoffich1ll995l) . SLIP does not assume that line 
scattering is depolarizing; it also does not rely on the 
Sobolev approximation, but instead performs full ra- 
diative transfer in regions of high optical depth. It is 
thus able to probe in detail the polarized line profiles 
that may arise from interaction with the circumstel- 
lar material. Another advantage to this method is 



that it can simulate emission not only from a central 
source but also from extended regions such as the 
warm CSM (see below). However, the code is still in 
the early stages of development, and does not include 
any Doppler effects from the expanding circumstel- 
lar material; this limits the extent to which we can 
compare model outputs to observed line profiles, but 
the stationary case is a useful first approximation, 
especially in cases of low CSM velocity. 

In the models presented here, a finite spherical 
source of photons at the center of the grid repre- 
sents the "photosphere" of the SN ejecta, while two 
scattering regions surrounding the ejecta represent 
the warm, stationary CSM and a shock-heated re- 
gion interior to the CSM, formed by its interaction 
with high- velocity SN material. Initially unpolarizcd 
photons are emitted from the surface of the photo- 
sphere with the synthetic Ha P Cygni profile from a 
Type IIP supernova, pro duced with the PHOE NIX 
stellar atmosphere code ( Hauschildt et al. 1999f ). In 
the hot shock region, I assume only a narrow line 
is emitted, with a line width of < 80 km/s; this is 
similar to the widths of the narrow hydrogen lines 
in the Type I l n SN 1997eg (7-40 km/s) observed 
by Salamanca ( 2002j ). I do not directly simulate 
heating of the CSM by emission from the supernova, 
but rather choose a CSM temperature and emit pho- 
tons from the volume of the region with the expected 
thermal continuum an d line spectra of hy drogen, as- 
suming Case B LTE ( Osterbrockl 119891) . All pho- 



tons then scatter within the CSM via electron scat- 
tering; the optical depth of the scattering region is 
also chosen as an input parameter. I also implement 
wavelength-dependen t free-free and free -bound ab- 
sorption effects as in IWood et al.l (11996 ). Pho tons 
within the Doppler core of the line (jLane Il999l ) ex- 
perience an additional bound-bound opacity; those 
absorbed by H atoms in this way are subsequently 
rc-cmittcd coherently and isotropically. 

When the code is run, SLIP emits photons se- 
quentially from the photosphere, shock region, and 
CSM and follows each until it becomes absorbed or 
exits the system; at each scattering event, the code 
updates the photon's Stokes parameters. Photons 
that exit the model system are binned by outgoing 
angle and the Stokes parameters are summed appro- 
priately in each bin. Thus each simulation produces 
a full three-dimensional model whose Ha flux and 
polarization spectra can be "viewed" from any de- 
sired direction. 

3. MODEL RESULTS 
I have created a grid of 144 models spanning two 
CSM geometries (ellipsoids and toroids of similar 
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Fig. 1. Representative sample of the model flux grid described in ij3] All models shown here have ellipsoidal circumstellar 
density distributions; none includes emission from the shock region. In the individual panel labels, the upper number 
indicates the CSM temperature in Kelvin and the lower indicates the CSM optical depth. Each panel spans a wavelength 
range from 5800 A to 7200 A; each spectrum is normalized to its own peak flux. Four profiles are plotted in each panel, 
representing viewing angles of 3°, 35°, 66°, and 89°. Model A (Table [1} is depicted in the middle right panel. 



sizes); CSM optical depths from 0.5 to 2; CSM lumi- 
nosities from 1-20% of the photospheric luminosities; 
and CSM temperatures from 10,000 K to 20,000 K. 
Emission from the shock region was either included 
at 10% of the photospheric luminosity or excluded 
completely. In order to keep computing times rea- 
sonable, these models included only 1.8 x 10 7 pho- 
tons each, enough to build up good signal in the flux 
spectrum but not the polarization spectrum. Since 
the goal was to match general features in both the 
flux and polarized flux, I used this grid to constrain 
parameter space for more computationally-intensive 
code runs including polarization. Figure[T]depicts re- 
sults for nine representative models in the flux grid; 
it shows that significant differences in the line pro- 
files arise for even small variations in optical depth 
and temperature of the scattering region. 

With the flux grid complete, I compared the sim- 
ulated Ha profiles with observed Ha profiles of Type 
Iln supernovae to determine which simulations to 
repeat for better signal in the polarized flux. This 
process is still underway, and the results will be pub- 
lished in an upcoming contribution. Here I present 
two representative simulations, Model A and Model 
B, that produce similar Ha flux spectra but have 
quite different polarization behavior. Each included 
1.6 x 10 9 photons, divided among 768 processors of 
the Seaborg parallel computing facility at the Na- 
tional Energy Research Scientific Computing Center 
(NERSC) at the Lawrence Berkeley Laboratory. Ta- 
ble [T] compares the parameters of these two models, 
while Figure [2] compares their Ha line profiles at a 



TABLE 1 

REPRESENTATIVE MODEL PARAMETERS 



Parameter 


Model A 


Model B 


CSM geometry 


ellipsoid 


toroid 


CSM optical depth 


1.0 


2.0 


Lcsm/ Lphot 


0.01 


0.1 


Lghock/ Lphot 


0.0 


0.1 


CSM temperature 


20,000 K 


15,000K 



range of viewing angles with the profile observed for 
SN 2000P (A. V. Filippenko 2004, private commu- 
nication) at day 13 post-discovery. Both models can 
reproduce the general observed line shape of a nar- 
row emission "spike" superposed on a broad base (I 
note that due to the SLIP code's limitation to sta- 
tionary scattering regions, the width of the broad 
line in these models arises solely from the input IIP 
spectrum arising from the model SN photosphere). 
While the profiles produced by Model A (the ellip- 
soid) are nearly completely degenerate in viewing an- 
gle, Model B (the toroid) shows a significant varia- 
tion with viewing angle, particularly in the strength 
of the "spike" relative to the broad base. Examina- 
tion of the other models in the grid suggests that the 
differences between these two families of Ha profiles 
are mainly due to the difference in CSM geometry 
between these two models. 

In Figure [3] I present the Ha line profiles of Mod- 
els A and B in polarized flux and compare them 
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Fig. 2. (top) De-redshifted Ha total flux profile of SN 
2000P at day 13 post-discovery (A. V. Filippenko 2004, 
priv. comm.). (middle) Simulated Ha profile arising 
from Model A (Table [TJ at varying inclination angles 
from the polar axis, (bottom) As in the middle panel, 
but for Model B. All line profiles have been normalized 
to 1 at the rest wavelength of 6563 A. 



Fig. 3. (top) De-redshifted Ha polari zed flux profil e 
of SN 1997eg at day 44 post-discovery (|Hoffmanll200r3 '). 
(middle) Simulated Ha polarized flux profile arising from 
Model A (Table H) at varying inclination angles from 
the polar axis, (bottom) As in the middle panel, but for 
Model B. All polarized line profiles have been normalized 
to 1 at the rest wavelength of 6563 A. 



with that of the Type Iln SN 1997eg (jLeonard et al 
2000bl : lHoffmar]|2006t ) at 44 days post-discovery. Re- 
call that polarized flux is percent polarization multi- 
plied by total flux; these profiles thus represent the 
spectra of the scattered light in each model. In po- 
larized light the Ha profiles look quite different than 
in direct light. The degeneracies that characterized 
Model A in direct light (Figure ^ have been lifted, 
raising the possibility of using polarized line profiles 
to diagnose inclination angle in cases of known ge- 
ometry. The two models now show line profiles quite 
distinct from each other; in particular, Model A pro- 



duces no narrow "spike" at the rest wavelength in the 
polarized flux, while Model B preserves the spike. 
Examination of the other models in the flux grid 
suggests that this difference is due not to the differ- 
ence in geometrical structure of the CSM between 
the two models, but rather to the presence of shock 
emission in Model B. Evidently the diffuse thermal 
emission from the CSM is sufficient to create a nar- 
row spike in direct light, but in polarized light this 
feature requires a more directional source of narrow- 
line photons (the shock region). Photons from this 
region are more likely to scatter and become polar- 
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ized while traversing the CSM than are photons that 
arise within its volume. However, Model B does not 
match the observed magnitude of polarization in SN 
1997eg, most likely due to its larger CSM optical 
depth. 

I continue to investigate these model results to 
pinpoint further diagnostics for geometry, tempera- 
ture, and optical depth of the circumstellar material 
and the presence of a narrow-line "shock" region in 
Type Iln cjecta. More detailed analysis will be pub- 
lished in an upcoming contribution. 

4. OBSERVED Ha LINE PROFILES OF SN 
1997EG 

As mentioned in my current SLIP code has 
the limitation of treating only stationary scattering 
regions, when in fact one expects velocity effects to 
be quite prominent contributors to line profiles in 
polarized light. Here I present two examples of line 
polarization effects in a Type Iln supernova that can- 
not yet be reproduced by the code but provide key 
information regarding the geometry of the circum- 
stellar material. Along with collaborators at San 
Diego State University and UC Berkeley, I have stud- 
ied the polarization spectrum of SN 1997eg (whose 
Ha polarized flux profile i s shown in Figu re El) ; our 
full analysis will appear in iHoffmar] ( 2007 ) . 

Figure [4] compares the Ha line profiles of SN 
1997eg in total flux (nearly all direct light) and scat- 
tered light at day 16, day 44, and day 93 post- 
discovery. Although the continuum polarized flux is 
only about 2% of the total at each epoch, I have nor- 
malized both spectra to the same scale for compar- 
ison of the line profiles. At all epochs the polarized 
lines are broader than the unpolarized lines. This 
implies first that the scattering region has a different 
geometry from that of the broad- line emission region, 
and second that the scattering region is expanding 
at a higher velocity than the emission region. The 
fact that the polarized lines are broader only in the 
blue wing suggests the redshifted side of the scatter- 
ing region may be self-occulted (or perhaps occulted 
by the SN ejecta) from our line of sight. If we postu- 
late a flattened toroidal or disk-like geometry for the 
scattering region, this result can help place limits on 
the spatial inclination of the CSM configuration. 

In Figure [5] I plot the polarized Ha lines of the 
three epochs in q-u space, a technique that allows 
visualization of all polarimctric information at once. 
Not only does the magnitude of polarization change 
across the Ha line at all epochs in SN 1997eg, but 
the position angle changes as well, in a manner that 
creates closed "loops" in the q-u plane (distinct 
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Fig. 4. Line profiles of the Ha line of SN 1997eg in to- 
tal flux (narrow smooth lines) and polarized flux (thick 
binned lines) at each of the three epochs of spectropo- 
larimetry (days 16, 44, and 93 post-discovery). Total 
flux spectra have been normalized to their respective line 
peaks. Each polarized flux spectrum has been binned to 
a resolution of 10 A and multiplied by the same normaliz- 
ing factor as its corresponding total flux spectrum, then 
by an additional factor shown in each frame to facilitate 
direct comparison of the line shapes. 

from "knots," which characterize a constant polar- 
ization signal with wavelength, and from straight 
lines, which arise from simple envelope expansion 
or line depolarization). This implies that the scat- 
tering region polarizing the Ha line has a different 
orientation than the one polarizing the continuum 
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Fig. 5. Polarization profiles of the Ha emission line of SN 
1997eg in the q-u plane for each of the three epochs of 
spectropolarimetry (days 16, 44, and 93 post-discovery). 
Data have been binned to a resolution of 50 A for clarity. 
Dashed lines represent negative (blueshifted) velocities; 
solid lines represent positive (redshifted) velocities. In 
each frame, the rest wavelength is shown with an open 
diamond and the origin of the plot is in the upper right 



light (presumably the SN ejecta). In particular, the 
"loop" shape implies that the symmetry axes of the 
two regions are different, and that the CSM occults 
the SN ejecta in such a way as to create an asymme- 
try in the Stoke s parameters across the line center. 
In the models of iKasen et al.l (|2003l ). such q-u loops 
are general features of two-axis systems; I postulate 



a toroidal geometry such as that shown in these au- 
thors' Figures 14 and 15. 

I note that the preceding results are both inde- 
pendent of interstellar polarization effects, for which 
the observed data have not been corrected. Simi- 
lar complex behavior in the H/3 and He I A5876 lines 
in the SN 1997eg spectrum is discussed fully in the 
upcoming article. The results combine to suggest 
that SN 1997eg is characterized by ellipsoidal ejecta 
that polarize the continuum light and a flattened, 
disk-like CSM exterior to the ejecta that polarizes 
the hydrogen lines. The key conclusion is that the 
ejecta and the CSM have different axes of symmetry. 

Multi-axis systems are becoming recognized as 
quite common in mass-loss scenar ios. P Cy 



(iNordsieck et al.ll200ll ; lMeaburnll200ll). n Car jSmit 
2006)), and VY CMa ([Humphreys! l2005f) are exair 
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pies of massive stellar systems in which the circum- 
stellar material shows evidence for multiple mass- 
loss episodes along different axes. Our observations 
of polarized line profiles in SN 1997eg suggest that 
the geometry of an asymmetric SN explosion may, 
in turn, be unrelated to the geometry of the progen- 
itor's stellar wind or its mass eruptions. Continued 
refinements to the SLIP radiative transfer code will 
allow me to construct more detailed models of the 
CSM surrounding SN 1997eg and related objects, as 
well as further probing the nature both of Type II 
SN explosions and the eruptions that characterize 
their progenitor stars. 
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entific Computing Center, US DOE Contract #DE- 
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Alex Filippcnko at UC Berkeley, Peter Nugent at 
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